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We study the ground-state properties of frustrated Heisenberg ferrimagnetic ladders with antifer- 
romagnetic exchange interactions and two types of alternating sublattice spins. In the limit of strong 
rung couplings, we show that the mixed spin-1/2 and spin-1 ladders can be systematically mapped 
onto a spin-1/2 Heisenberg model with additional next-nearest-neighbor exchanges. The system is 
either in a ferrimagnetic state or in a critical spin-liquid state depending on the competition between 
the spin exchanges along the legs and the diagonal exchanges. 
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I. INTRODUCTION 

In the past years, there has been increasing the- 
oretical interest in the quantum ferrimagnetic sys- 
tems since the experimental realization of bimetallic 
quasi-one-dimensional (quasi-lD) magnetsii^i^. Pi- 
oneering experiments on synthesizing the bimetaUic 
chain compounds with each unit cell containing two 
spins were carried out successfully by Kahn et. ali. 
Typical compounds include two families of ferrimag- 
netic chains described by ACu{pha){H20)z-nH20 
and ACu{pbaOH){H20)3.nH20, where pba=l,3- 
propylenebis(Oxamato), pfeaOi7=2-hydroxo-l,3- 
propylenebis(Oxamato) and A=Ni, Fe, Co, and 
Mit^i^. These materials are quasi-lD bimetallic molec- 
ular magnets containing two different transition-metal 
ions per unit cell alternatingly distributed on a chain. 
Most of these materials are described by Heisenberg 
mixed-spin models with antiferromagnetic interactions. 
This has stimulated theoretical studies on the mixed-spin 
systems. A number of recent studies have been focused 
on ground-state properties of mixed-spin chain^'^i^'^ 
and ladders^ i^°i^^i^^i^'^i^'^i^^ composed of two kinds of 
spins with half-integer and integer spins alternatively. It 
is well-known that a half-integer antiferromagnetic spin 
chain has gapless excitation spectrum and the integer 
ones show the Haldane gap in their low energy spectrum, 
however the ID ferrimagnets behave differently and 
exhibit intriguing quantum spin phases and thermo- 
dynamic properties. One of the intriguing features of 
the mixed-spin chain is that the ground state has a 
long-range ferrimagnetic order. 

While the physics of the mixed-spin chain is well un- 
derstood, the spin ladders are still under intensive inves- 
tigation since they exhibit rich phase diagrams. So far, 
a large class of mixed-spin two-leg ladders are investi- 
gated, but most of studies concern the physics of unfrus- 
trated mixed-spin ladders. For the unfrustrated bipartite 
ladder with alternating half-integer and integer spins on 
neighboring sites, the ground state is shown to exhibit 
long-range order—. In general, for the uniform spin lad- 
ders, frustration reduces the antiferromagnetic correla- 
tions and may produce various exotic quantum ground 
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FIG. 1: Schematical representation of the mixed-spin ladder 
with diagonal exchanges. The open and close circles represent 
spins s. —1/2 and 5^,2 ~ 1 respectively. 



states such as the dimerized statei^iiiii^. In comparison 
with the uniform systems, less attention is paid to the 
role of frustration in the mixed spin systems and it is 
quite interesting to explore how the ferrimagnetic order 
is affected by the frustratio n^^^d^ ^ Recently, the effect of 
magnetic frustrations due to diagonal exchange bonds in 
a mixed-spin ladder was studied by using the approach of 
exact numerical diagonalization^^ for a special case with 
equal strengthes of the rung coupling and diagonal cou- 
plings. It was shown that the long-range ferrimagnetic 
state may be destroyed and a singlet state appears at a 
large value of the frustration parameter. However, it is 
not clear how the diagonal frustrations compete with the 
exchanges along the legs and whether an exotic interme- 
diate phase exists for arbitrary rung coupling. Further- 
more, an analytical analysis which works in both regimes 
with weak and strong frustrations is still lack. The gen- 
eral analytical method based on the spin wave theory is 
not capable to deal with the regime with strong frustra- 
tion where the long-range order is destroyed. In order 
to provide a general understanding to frustration effect 
in the ferrimagnetic ladder system, in this work we carry 
out analytical investigation on the extended antiferro- 
magnetic mixed-spin ladders in the limit of strong rung 
couplings. When the rung coupling is dominated, it is 
convenient and natural to map the original model onto 
an effective spin-1/2 exchange model. Based on the ef- 
fective model, we show that the ground-state properties 
for both the frustrated and non-frustrated ladders can 
be understood within the same framework and the phase 
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transition from a ferrimagnetic state to a critical spin- 
liquid state is also discussed. In our theory, we treat 
the weak and strong frustrations on same footing. For 
generality, we consider the most general frustrated lad- 
der model with arbitrary strengthes of exchanges of spins 
along the two legs as well as the diagonal exchanges. Our 
model can interpolate between a variety of systems, ex- 
hibiting remarkably rich ground state behavior with both 
ordered and disordered phases. 



II. MODEL AND THE EFFECTIVE 
HAMILTONIAN 

As shown in Fig.l, the extended frustrated mixed-spin 
ladder is made of two types of spins with magnitude 
s.^ = \/2 and 5*^,2 = 1 located on the rungs of the lad- 
der. The corresponding Hamiltonian of the mixed-spin 
ladder takes the form 



H = Ho + H' 



(1) 



the one with Jx = J'x = is equivalent to the railroad 
ladder model; For Ji — J2 and Jx = J'x^ our model re- 
duces to the ladder model investigated in Refii^, where 
only the case with J± = Jx was considered; For Ji = J2 
and Jx=0, the model is the dimerized zigzag mixed-spin 
ladder^iii. 

In this work, we will focus on the strong coupling limit 
with Jj_ ^ Ji, J2, Jx, Jx- In this limit, the interac- 
tions between the neighboring rungs can be treated as 
perturbations of the system of uncoupled rungs. It is in- 
structive to start by considering the two-site problem on 
an isolated rung with the Hamiltonian given by 

Hi = J±Si^i ■ Si^2- 

Eigenstates of the local Hamiltonian on a rung can be 
classified according to the value of the total rung spin Si. 
The two spins with magnitude Si^i = 1/2 and Si^2 = 1 
can combine into Si — 1/2 and 3/2. It is easy to get 
the eigenenergy E1/2 = — J± and E^/2 J±/2. The 
corresponding eigenstates are 
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Hi = j± ^ • Si^2 
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where 
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J2<5'.i,2 ■ 'S'i+1,2 
JxSi,l • '5'i+l,2 
J'y, 5*1,2 • + 



(2) 



(3) 



(4) 



(5) 



Here N is the number of rungs, and Si,i and 5^,2 repre- 
sent the spin-1/2 and spin-1 operators respectively. For 
the antiferromagnetic ladder, all the coupling parameters 
J±, Ji, J2, Jx, J'x > 0- The intrachain couplings along 
the up and down legs are denoted by Ji and J2 respec- 
tively. The interchain coupling across the rungs is J± and 
the diagonal exchanges between the rungs are Jx and J'^ . 
The summation is over the length of the chains and the 
periodic boundary is assumed. For such a ladder model, 
the frustration comes from the competition between the 
intrachain couplings (Ji and J2) and the diagonal cou- 
plings (Jx and J'y). The model of ^ covers a variety 
of known models. The model with Ji = J2 = corre- 
sponds to a bipartite mixed-spin ladder model whereas 



D^) = ^ 



1 

1 

71 



V2 



-\.0) -V2 



1 

1 



for the doublet and 
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for the quartet, where \sz,Sz)i = \sz)i 1 ® \Sz)i 2- It is 
obvious that the ground state of a rung is a doublet and 
the excited state is a quartet with an excitation energy 
gap 3J_l/2. Therefore in the strong rung-coupling limit 
spins on each rung of the ladder favor forming a doublet. 
Since each rung can be either in the doubly degenerate 



state 



Di 



■ ^-i / I the ground state of the zero-order 

Hamiltonian Hq is 2^— foldly degenerate. When we go 
beyond the two-site problem, we need consider the inter- 
rung couplings. In general, the inter-rung exchanges will 
lift the degeneracy of the zero-order ground state and 
leads to an effective Hamiltonian acting on the ground- 
state Hilbert space of Hq. 

We then derive the effective Hamiltonian of the original 
ladder model in the truncated Hilbert space composed of 
product of rung doublets by using perturbation method. 
Similar schemes have been applied to study uniform spin 
ladders in the strong-coupling limit^^^. To the first or- 
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der and up to a constant of —NJ±, the effective Hamil- 
tonian can be represented as 

^e// ^ X! Hi,i+1 \Vi,i+l) {ViJ+l\ (6) 

i 

where , = 



D 



I Da_i ) I are four- 

foldly degenerate. It is convenient to introduce pseudo- 
spin-1/2 operators % which act on the states and 
1-0-1/2)^ and are defined as 



D_ 



1 

±2 

D_ 



D_ 



0. 



With the above notation, we can identify the following 
relations 
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Therefore we can rewrite terms of | in terms 

of the pseudo-spin operators. For example, we have 



for = 



D_ 



and 



2 

given by 



i+l 



and the corresponding coefficient is 



2 18 

After some algebras, we can rewrite the effective Hamil- 
tonian ([6]) as the following form 



(1) 



1 



i 'i+l 



H 'i+l 



with the effective parameter given by 



(7) 



J. 



(1) 



4/Jx + (8) 

9\J+ J+J 9J+ 9 J_L ■ ^ ' 



Next we carry out calculation of the second order per- 
turbation which is expected to give correction to the ef- 
fective coupling parameter and the ground energy. In the 
first order calculation, the higher excited states (quar- 
tets) on the isolated rungs do not play role. However, 
such rung quartets give contribution to the higher order 



correction. The second order correction can be described 
by 

^ I Hij+i \m) (m| Hi^i+i \vi^i+i) 



H 



(2) _ 



i,7n^Q 



En ~ E„ 



(9) 



where Eq — 2Ei/2 = — 2J^ and |m) are the inter- 
mediate states with at least one quartet residing in 
the neighboring i th and (i -I- l)th rungs. Therefore 

E,n = + = -J±/2 if |m) = {\Da), |Q/3),+i 

or |Q/3)J-Da)i+i}, and E,n = 2£^3/2 J± if |m) = 
IQ/3)JQ/3'),+i where a = -1/2,1/2, and /3,/3' = 

-3/2,-1/2,1/2,3/2. After tedious but straightforward 

(2) 

calculation, we can represent H^^j in terms of the 
pseudo-spin operators 



Hfff/J+ 



J2 [jefA ■ n+1 + c' 



(10) 
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and 
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The second order correction can also produce a term due 
to the three-site process which is described by 

(M»,i-|-l,i+2| -ffe// \'^i,i+l,i+2) 

{fJ-i,i+l,i+2 \ Hi^i+i \m) (m| HiJ^i,i+2 \uiA+l,i+2) 
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£^0 — Em 

{^J'iA+l,i+2\ Eli+l,i+2 [m) {m\ iJi,i+l \vi^i+l.i+2) 
En — Em. 



where Eq — 3£'i/2 = — 3Jj^ and \m) are the intermedi- 
ate states with a quartet residing in the {i -\- l)th rung. 
In terms of the pseudo-spin operators, the Hamiltonian 

H^jf is finally simplified to 
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(14) 



Therefore, up to the second order, the effective Hamil- 
tonian of the original ladder can be written as 



N 



Heffl JA-=^ {JeffTi ■ Ti+i + J'ef fTi ' n+2 



(15) 



1=1 
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where J^ff = J^^j + J).jf and c = — 1 + c'. The effective 
Hamiltonian (jlSp describes a spin- 1/2 Heisenberg chain 
with additional next-nearest-neighbor (NNN) exchanges. 



III. PHASE DIAGRAM AND DISCUSSION 

It is convenient to determine the phase diagram of 
the original ladder model from the effective Hamiltonian 
(|15p since it has a much simpler form and was widely 
studied. Since the terms with NNN exchanges come 
from the second-order perturbation, in general we always 



have 



|Je//|. A small NNN ferromagnetic in- 
teraction is an irrelevant perturbation to the spin chain 
modeliii^. Therefore we may approximate to determine 
the phase diagram by omitting the term of J^y ^ . The ef- 
fective Hamiltonian ([TS]) with J'^^j — describes either 
a ferromagnetic or an antiferromagnetic spin-1/2 Heisen- 
berg chain depending on the sign of the effective param- 
eter. Based on the effective Hamiltonian, one is ready to 
get the ground energy of the mixed-spin ladder ([T]) in the 
strong coupling limit: 



Eg/NJ^ 



--H 



Se' 



c, 



(16) 



where eh is the ground energy per site of an isotropic 
spin-1/2 Heisenberg chain and Se' is the energy correction 
due to the NNN exchanges which reads 



6e' = {GS\H'}J] \GS) /NJa^ 

with \GS) representing the ground state of Heisenberg 
chain. From the corresponding results of the well-known 
Heisenberg chai n^^'^"^ , we have 



£H/Jeff - 



1/4, 
1/4- 



ln2, 



for Jeff < 
for Jeff > 



When the system is in the ferromagnetic state, it is ready 
to obtain 

Se' ^ J'eff/A. 

When the system is in the antiferromagnetic phase, no 
an analytical result is available. 
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FIG. 2: Phase diagram of the mixed-spin ladder with Ji = J2 
and Jx = J'x in the strong coupling limit. The solid line 
is the phase boundary determined by Je// = whereas the 
separated dots denoted by stars are obtained from the exact 
diagonalization of the spin ladder model. 



From eq.®, we see that the diagonal exchanges tend 
to produce an effective ferromagnetic coupling whereas 
the exchanges along the legs lead to an effective antifer- 
romagnetic coupling. Both the ferromagnetic and antifer- 
romagnetic spin-1/2 Heisenberg chains are exactly solved 
by the Bethe-ansatz method^i^, but they have very 
different ground-state properties. The ferromagnetic 
ground state consists of all spins parallel and thus has 
long-range order, whereas the antiferromagnetic ground 
state is a spin singlet which can be described by the crit- 
ical spin-liquid phase. If Je// < 0, the effective model is 
a ferromagnetic chain which exhibits gapless excitations 
and ferromagnetic ground state. Equivalently, within the 
first order approximation the original mixed-spin ladder 
has ferrimagnetic ground state if Jx + J'^ > Ji/4 -I- 4J2. 
On the other hand, the effective model is an antiferro- 
magnetic chain if Je// > 0. Correspondingly, the ground 
state of the original mixed-spin ladder is a critical spin 
liquid if J^ + J'^ < Ji/4 -I- 4J2. 

As we have mentioned that our mixed-spin ladder 
model includes a series of models as its special cases, 
we may understand the ground-state properties of these 
models by using our derived effective Hamiltonian (jlSp . 
When Ji — J2 = 0, our model reduces to a bipar- 
tite mixed-spin ladder model. The effective model 
is a ferromagnetic Heisenberg chain because we have 
Jeff < 0. Therefore we can conclude that the total spin 
of the ground state is N/2. When Jx — J'^ ~ 0, we get 
the railroad ladder model composed of coupled spin-1/2 
and spin-1 chains. It is ready to conclude that its ground 
state is a singlet because the effective model ([T5| is an 
antiferromagnetic Heisenberg chain with Je// > 0. For 
the frustrated spin ladder with Ji = J2 and Jx = J'x, 
it is easy to see from the effective model (IT5|) that there 
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exists a quantum phase transition arising from the com- 
petition of Ji and Jx • The transition point is determined 
by Jef / = with 

8Jx^l7Ji 56 /Ji JxV 

From the solution of the above equation, we can get the 
phase diagram of the mixed spin ladder as shown in Fig. 2. 
Above the phase boundary the original ladder is in a 
fcrrimagnetic phase corresponding to Jeff < 0, whereas 
the system is in an antiferromagnetic phase with Jef / > 
below the phase boundary. 

In order to check how good is our theory based on 
the perturbation method and the validation of the phase 
boundary determined by Jeff = 0, we study the origi- 
nal spin ladder model by exact diagonalization method 
and determine its phase diagram numerically. The phase 
boundary determined by the effective spin chain model 
is not sensitive to the size of the system (Actually the 
phase boundary is completely independent of the size 
when we omit the term of NNN interaction) . Therefore 
for the purpose of verification of our analytical result, 
it is enough for us to consider a 2 x 4-size ladder which 
can be diagonalized by the exact diagonalization method. 
The phase boundary can be determined simply by the 
ground state degeneracy of the system. In comparison 
with phase boundary obtained by omitting the term of 
the NNN exchanges, we find that they agree very well in 
the strong coupling limit and as expected, the analytical 
result begins to deviate the exact numerical result when 
Ji/J± increases. 

Finally, we address a special point with Ji = J2 = 
Jx = J'x where the total spin on a rung is a good quan- 
tum number. From eq.®, it is obvious that the model is 
effectively described by an antiferromagnetic spin chain 

with an effective coupling constant of Ji and the high- 

(2) 

order corrections vanish due to J^A — J'^j^ = 0. Actu- 
ally this is an exact conclusion for this special case. This 
is clear if we rewrite the original Hamiltonian as 



N 



Ji 



Jl^i ■ Si+i 



IIJ, 



where = S,iS^ + 1) with S, = 1/2 or 3/2 and S, 
is defined as Si = Si,i -I- 5*^,2 denoting the total spin on 
the ith rung. In the strong coupling limit J± 3> Ji, J_l 
forces Si to take the value of 1/2 and thus the ground- 
state properties of H is described by an effective spin-1/2 
antiferromagnetic Heisenberg chain. From eas. pTj) and 



(fT^ . we see that J^ 



(2) 
eff 



when Ji — J2 — Jx 



J'x ■ This implies that the higher order correction also 
conforms to the exact result in this special case. 

IV. CONCLUSIONS 



In conclusion, we have studied the ground-state prop- 
erties of a generalized mixed-spin ladder with diagonal 
exchanges in the limit of strong rung couplings. By map- 
ping it to an effective spin-1/2 Heisenberg chain with 
additional NNN exchanges, we find that the diagonal 
exchanges lead to the ferromagnetic effective coupling 
whereas the exchanges along the legs produces the an- 
tiferromagnetic effective coupling. The ground state of 
the effective Hamiltonian is either ferromagnetic or an- 
tiferromagnetic depending on the competition between 
these two opposite processes. With the help of the effec- 
tive Hamiltonian and omitting the additional small NNN 
exchanges, it is straightforward to analytically determine 
the transition point from the fcrrimagnetic phase to the 
critical spin liquid phase. The phase boundary is found 
to agree with the result obtained by exact numerical diag- 
onaliztion of the original spin ladder model. Our results 
show that the strong coupling approach provides a sim- 
ple and unifying way to exhibit the rich physics of the 
mixed-spin ladders. 
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